DESY 08-151 
Edinburgh 2008/39 



Non-Perturbative Renormalization of Three-Quark 

Operators 

Meinulf Gockeler ^, Roger Horsley ^, Thomas Kaltenbrunner ^, 
Yoshifumi Nakamura ^, Dirk Pleiter ^, Paul E. L. Rakow ^, 
Andreas Schafer ^, Gerrit Schierholz Hinnerk Stiiben ^, 
Nikolaus Warkentin ^, and James M. Zanotti ^ 

(QCDSF/UKQCD Collaborations) 

^Institut fiir Theoretische Physik, Universitdt Regensburg, 93040 Regensburg, Germany 

^School of Physics and Astronomy, University of Edinburgh, Edinburgh EH9 3 JZ, UK 

^Deutsches Elektronen-Synchrotron DESY and John von Neumann Institutfur Computing 

NIC, 15738 Zeuthen, Germany 

'^Theoretical Physics Division, Department of Mathematical Sciences, University of 
Liverpool, Liverpool L69 3BX, UK 

^Konrad-Zuse-Zentrum fiir Informationstechnik Berlin, 14195 Berlin, Germany 



Abstract 

High luminosity accelerators have greatly increased the interest in semi-exclusive and ex- 
clusive reactions involving nucleons. The relevant theoretical information is contained in 
the nucleon wavefunction and can be parametrized by moments of the nucleon distribu- 
tion amplitudes, which in turn are linked to matrix elements of local three-quark operators. 
These can be calculated from first principles in lattice QCD. Defining an RI-MOM renor- 
malization scheme, we renormalize three-quark operators corresponding to low moments 
non-perturbatively and take special care of the operator mixing. After performing a scheme 
matching and a conversion of the renormalization scale we quote our final results in the MS 
scheme at /i = 2 GeV. 
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1 Introduction 



Distribution amplitudes play an essential role in the investigation of the internal nu- 
clear structure. Exclusive high-energy processes can be factorized into hard and soft 
subprocesses, where the hard subprocess can be evaluated perturbatively and is char- 
acteristic for the reaction in question. The soft subprocess, described by the nucleon 
distribution amplitude contains the information about the distribution of the 

three valence quark momentum fractions Xi inside the nucleon [|1I2I3I4I5| . 

The great interest in this quantity stems from its importance for, e.g., the calculation 
of the electromagnetic form factors of the nucleon and their scaling behavior. These 
form factors describe a nucleon absorbing a virtual photon of squared momentum — 
while remaining intact. According to [2], for Q"^ oo, the magnetic form factor of 
the nucleon can be written as a convolution of three amplitudes: first, the distribution 
amplitude $ for finding the nucleon in the valence state with the three quarks having 
definite momentum fractions Xi, second, the hard scattering kernel Th, which describes 
one of the three quarks absorbing the photon, and finally the complex conjugate of $ 
that gives the amplitude for the outgoing quarks to form a nucleon again: 

G'M(g') = Jjdx] Jjdy]^*{y„Qy)TH{x,,y,,QMx,,Q^) (l + 0{m%/Q^)) . 

(1) 

Here [dx] = dxidx2dxsS{l — Y^i^i), Qx = niinj(xj(5) and denotes the nucleon 
mass. 

Distribution amplitudes are genuinely non-perturbative quantities. Hence they are in- 
accessible for perturbation theory and must be calculated by other means, e.g., by 
lattice QCD. After a pioneering study in the late 1980s [6] only recently the first quan- 
titative results from lattice QCD have been published [7 8 9|. After performing an 
expansion near the lightcone, moments of the nucleon distribution amplitudes are ex- 
pressed in terms of matrix elements of local three-quark operators that are evaluated 
between a nucleon state and the vacuum. These three-quark operators typically con- 
sist of a combination of covariant derivatives acting on the three quark fields, which 
are located at a common space-time coordinate x. Furthermore, the operators for the 
nucleon are color singlets and have isospin 1/2. 

As the three-quark operators pick up radiative corrections and are subject to mixing 
with other operators, their renormalization is a vital ingredient for any lattice calcu- 
lation. In a previous paper [fTOll we have derived irreducible multiplets of three-quark 
operators with respect to the spinorial hypercubic group H(4) and have discussed the 
mixing properties of these operators in detail. The present work will focus on the non- 
perturbative renormalization of these three-quark operators in the isospin-1/2 sector. 
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The paper is organized as follows. In the first section, we will introduce an RI-MOM 
renormalization scheme which will be implemented in our lattice calculations. Then 
we will explain how to perform a perturbative scheme matching to MS and derive the 
anomalous dimensions of the operators in question. The following two sections will 
focus on a discussion of our results for the renormalization matrices, which are finally 
quoted at /i = 2 GeV in the MS scheme. In the final section we will demonstrate how 
to renormalize moments of the nucleon distribution amplitude and present the results 
of consistency checks between the renormalized moments. 



2 Lattice Renormalization 

Our main aim is the derivation of non-perturbative renormalization coefficients for 
three-quark operators in lattice QCD and their application to the renormalization of 
moments of the nucleon distribution amplitude. As discussed in the introduction, a 
subsequent convolution with the hard scattering kernel leads then, e.g., to estimates 
for the electromagnetic form factors of the nucleon. Therefore the renormalization of 
the hard scattering kernel and the three-quark operators inside the nucleon distribu- 
tion amplitudes must be carried out in a consistent way, i.e., the same renormalization 
scheme has to be applied. As the perturbative results for the hard subprocess are usu- 
ally given in the MS scheme, the distribution amplitudes must also be renormalized in 
this scheme. 

On the lattice, however, it is not possible to implement the MS scheme. Therefore 
one introduces the RI-MOM renormalization scheme [1 1 J. It is applicable both on the 
lattice and in the continuum. Then one first renormalizes the matrix elements on the 
lattice within this scheme. Afterwards one applies continuum perturbation theory to 
calculate a matching function between both schemes and extracts non-perturbatively 
renormalized lattice operators in the MS scheme. In this section we will discuss the 
setup and implementation of the lattice renormalization. 



2.1 The Three-Quark Operators 

We are interested in three-quark operators for the nucleon. Choosing the flavors u, u 
and d for definiteness, their general form is 

■ {D,, . . . D^Mx)^,)^, {D,, . . . D^^d{x)y)^, , (2) 
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with analogous definitions for Of'^'^{x) and Of'^^{x). After projection onto isospin 1/2 
one ends up with operators for the proton, and subsequent exchange of the u and d 
fields leads to neutron operators. The color indices c'j will be suppressed when not 
required. The coefficient tensor T*^*^ determines the symmetry structure of the spinor 
indices a', (3', 7' and of the space-time indices fXj, uj, Xj of the co variant derivatives D 
contributing to the operator Oj. 

In [[TOB we have derived multiplets of three-quark operators that transform irreducibly 
under the spinorial symmetry group of the hypercubic lattice H(4). These operators 
reduce the problem of mixing under renormalization to mixing among multiplets be- 
longing to equivalent representations of H(4). Hence they define our choice for the 
coefficient tensors T^^\ Appendix B of [lOl contains all linearly independent sets of 
potentially mixing isospin- 1/2 three-quark operators with leading twist and up to two 
derivatives. They provide the basis for the (non-perturbative) renormalization. Opera- 
tors with total derivatives are automatically taken into account due to the identity 
dfiifgh) = {D^f)gh + f{D^g)h + fg{D^h) for any color-singlet operator made of 
the three quark fields /, g and h. Note however that this continuum relation holds only 
up to discretization errors on the lattice. 

Let us finally give an example of a typical three-quark operator belonging to an irre- 
ducibly transforming multiplet. We have, e.g., the following operator with two deriva- 
tives: 

- {Da)^\y{Da)\yd\''^ + {Da)^\^d\Da)\yu\''^ 

- 2 ■ {Da)^\y{Da)\yd\''^ +2 ■ {Da)^\^d\Da)\yu\''^) . 

(3) 

The superscript MA (standing for "mixed antisymmetric") indicates that the opera- 
tor has isospin 1/2 while the subscript fg means that the derivatives act on the first 
and second quark, compare [10]. Pauli matrices are used to contract the co variant 
derivatives. The curly braces denote independent total symmetrization in the dotted 
and undotted indices of the Weyl representation. 

2.2 Calculational Method 

We introduce a correlation function for the non-perturbative renormalization of the 
three-quark operators by contracting the operators with three external quark sources, 
namely u{zi), u{z2) and ^(,23): 

I(a{zi)acji{z2)fic^{z^)^c30i{x)). (4) 
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We have two w-quark lines and one d-line running from the three space-time coordi- 
nates Zi into a vertex at x. There they are connected according to the coefficient tensor 

We proceed in analogy to the case of quark-antiquark operators II11I12II . We impose 
fixed momentum on the external quark lines and evaluate the correlation function in 
momentum space. The result is a four-point function G{pi,p2,P3)l^j^^ that depends on 
the momentum of the three external quark lines and carries three spinor indices: 

G{pi,P2, P3)%^ = ^ exp(+ipi ■ Zi + ip2 ■ Z2 + ip3 ■ Zs) 

■ exp(-i(pi +P2+ P3) ■ x) eciC2C3 {^zi)acM^2)f3c2d{z3)yc30i{x)). (5) 

This four-point function is obtained as the ensemble average of quark field contractions 
on the individual gauge configurations. Let us denote such a contraction on a single 
configuration by the brackets [...]. For Oi = O^^'^ we then have, e.g.: 

G{pi,P2,P3)\^(i^ = X! ^c;44 eciC2C3 ^/li^.../i^i/i...i/„Ai...A,a'/3'7' 
config. 



■ — ^ exp(-z(pi +P2+ Ps) ■ x) exp(+ipi ■ zi + ip2 ■ Z2 + ips ■ Z3) 

^ X,Zl,Z2,Z3 




where is the number of gauge field configurations. Analogous expressions hold for 
the four-point functions of the operators (9""^" and Of"^^ as well as for the operators 
with definite isospin. 

Let us now define 

K{x,p)aj^ci,a2C2 = Yexp{ip ■ z)[u{ 

]• (7) 

z 

This quantity can be determined on every single gauge configuration by inverting the 
massive Dirac operator M of the action on a momentum source [lT2l|: 

co,aiCi {y,x)K{x,p) = exp{ip ■ y)5aoa2^coc2- (8) 

Rewriting the above four-point function G in terms of K provides insight into the 
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lattice implementation: 

config. 

+ K • • • D;^),^,, ^(^5Pi);'4',«c. ■ K • • • D:^,,, K{x,p,)l,,,^,^^ . (9) 

The most expensive step in the calculation is the evaluation and symmetrization of 
the spin-color combinations which is, naively speaking, due to the existence of an 
additional quark and antiquark field with 12 spin-color indices, by a factor of 12 x 12 = 
144 more expensive than for mesonic operators. So, special care has to be taken when 
implementing these contractions. As the correlation function is not gauge invariant, all 
configurations are gauge fixed to Landau gauge. 

In a final step we amputate the external quark lines of the four-point function G^*) to 
arrive at the three-quark vertex F*^*^: 

^(^1,^2,^3)^7 = r(Pi 

) P2,Pz)\]l3i-j'S{j)i)a'cxS{p2)/3'(3S{j)'i).yi^, (10) 

where the quark propagators are defined by 

S{p)aia2 = ^I]("(^)"i'"(?/)"2)exp(-ij9- (a;-?/)). (11) 

Note that we can reuse the quantity K{x, p) introduced in (|7]) to calculate the propaga- 
tor ^(p): 

1 1 

S{p)a^a2 = Yl I]^(a;,p)aia2exp(-ip-x). (12) 

''^ config. X 

In the following section we will define our renormalization scheme based on the vertex 



2.3 An RI-MOM Renormalization Scheme 

We want to introduce a renormalization scheme that is applicable on the lattice and 
in the continuum. For quark-antiquark operators such a scheme has been proposed in 
[fm and is widely known as the RI-MOM scheme. In [il3ll it was also used for the 
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renormalization of proton decay matrix elements. To study the mixing of our three- 
quark operators with up to two derivatives, we have slightly modified this approach. 
In the following, we will set up our modified RI-MOM renormalization scheme. As, 
in general, mixing is a central issue of the renormalization, we will consider mixing 
matrices explicitly from the very beginning. 

The renormalized counterpart of a general regularized operator Oi is given by 

or = z,p,. (13) 

Here Zi^ is the renormalization matrix and operator mixing shows up in non- vanishing 
off-diagonal elements of Z. 

We will define the renormalization matrix Z™^^ for three-quark operators by projec- 
tions of the lattice-regularized three-quark vertex F introduced in eq. (flOl) . In the fol- 
lowing we will distinguish between the tree-level, the lattice regularized (lattice spac- 
ing a) and the renormalized vertices rf^*^, F'^" and r™"^, respectively. Let us now in- 
troduce a set of projectors Pk in spinor space that fulfill the following orthogonality 
condition with the tree-level vertices: 

P^VT[pX.V2,V■^=h^. (14) 

At some renormalization scale fixed by the mean squares of the three external quark 
momenta we then require the renormalized three-quark vertex to fulfill the same equa- 
tion. This yields the renormalization condition 

Pfcrf^^(pi,P2,P3;/i)l^2=^^p2/3 = bki- (15) 

With r^'*^ = z5'"'^Tf ' we can introduce the auxiliary variable Z^- for the renor- 
malization of the vertex: 

{Zvlm)^A^) = Pj^T\Pi,P2,P3)\,^=j:^p^/3 (16) 

In any scheme Z^ is related to the renormalization matrix of the three-quark operators 
O by the quark field renormalization Zg. To compensate for the amputated quark legs 
one needs a factor of Z^"^ for each of them: 

z,,(^) = ZMY^' z,iP2)'^' z,ip,Y/' . 4(/i). (17) 

This fixes the three-quark operator renormalization matrix in the RI-MOM scheme: 
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As usual we determine the factors Zq in the RI' scheme [fTTl . The peculiarity of the 
mRI scheme lies in the special definition of the projectors for the three-quark vertex F 
and will be explained in the next subsection. 



2.4 Choice of the Projectors 



In order to determine the RI-MOM renormalization matrix Z we have introduced but 
not yet defined a set of projectors Pj. The only restriction is given by eq. (fT4l) . 

We now turn to vector notation. Fj is a tensor of rank three in spinor-space and can be 
interpreted as a vector f of dimension 4^. Then we can interpret the projectors as 
orthogonal projections onto vectors vp^ of the same dimension: 

PfcF,= (t;p„t;r,), (19) 

where 

(i^i,^2) = E(^i)r (^2),- (20) 
3 



The task is now to construct a set of vectors vp^ that fulfills the normalization condition 
eq. (fT4l) . which reads in vector notation 

(fp^jWrtree) = (21) 

We choose the vectors vp^ as follows. We start with an auxiliary vector 

v'p^ = urfs (22) 

and project it onto the orthogonal complement of the space spanned by the vectors 
Wptree, j ^ k. Thls rcsults in an altered vector Vp^. Taking care of the normalization we 
finally define 



1 ,, 



By now all constituents of the mRI renormalization scheme are defined. We summarize 
the method by rewriting eq. (fTSi) in vector notation: 



(24) 

Note that the above renormalization condition will in general depend on the geometry 
of the external momenta, i.e., the angles between the four-momenta p^. In the end this 
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dependence will be cancelled by the scheme matching. We want to stress furthermore 
that, due to its general structure, the method is not limited to the case of three-quark op- 
erators and four-point functions discussed in this paper, but is applicable to the general 
class of n-point functions. 



3 Perturbative Calculations 



In the introduction we have emphasized the importance of renormalized three-quark 
operators for nucleon distribution amplitudes. To calculate observables, both the distri- 
bution amplitude and the hard scattering kernel must be given in the same renormaliza- 
tion scheme. In the following we will explain, how a matching of the RI-MOM scheme 
to MS can be achieved with the help of continuum perturbation theory. Moreover we 
will study the dependence of the renormalization coefficients on the renormalization 
scale fi. As for the lattice computation, also here all calculations have to be carried out 
in the Landau gauge. 



3. 1 Scheme Matching to MS 



Let us start by writing eq. (fT3l) explicitly in both renormalization schemes: 

Of^ = zfo,, (25) 

with O™^^ and Of^^ denoting the renormalized operators in the modified RI and the 
MS scheme, respectively. If we introduce the scheme matching matrix 



^MS^mRI^^MS(^(^mRI)-lj^^^ (26) 

we get a relation between the operators renormalized in the MS and mRI schemes: 

^MS ^MS^mRI^mRI (27) 

In the following we derive the matching functions with dimensional regularization in 
continuum perturbation theory. 

To this end we proceed with a one-loop perturbative expansion of Z^^ and Z™^^. Up 
to (9(e) the renormalization matrix reads in the mRI scheme: 

Z^^ = 6,, + ^ (Z^''') • (28) 



9 



Here we have used the renormalized strong coupling ^^(/i) = gR{jj,y /An and adopted 
the following conventions for the dimensional regularization: 

111 1 

- =- + -log47r--7B, (29) 
e = 4 - ci. (30) 

In the numerical evaluation of as we use A*^^ = 261 MeV [fT4l . Comparing eq. (|28] ) 
with the analogous expression for Z^^ and noting that the conversion between both 
must be finite, the scheme matching matrix in first order becomes 

Zf--^if^) = 6., - ^ {zr),^^ + Oial). (31) 



This means that 2^!^s^mRi derived from the perturbative expansion of Zl 



alone. In the following section we will discuss this in more detail. 



MS^mRI 



3.2 Determination of Z 



According to eq. (fTTI) the renormalization of the three-quark operators consists of four 
parts: 

zf' = {zfy/' {zfy/' (zf ■ zi^^^^. od 

The scheme matching can be performed independently for each renormalization factor 
so that 

In two-loop order the matching of the quark field renormalization reads IfTSl : 



4^f ^-^ + 12C3^-26e+-., + 12C3--j, 

(34) 

where ^ is the covariant gauge parameter, compare Appendix |Al The scheme matching 
for the renormalization coefficient Z^ of the three-quark vertex will be determined in 
one-loop order. In analogy to eq. (|3TI) we have 



,7r,MS^mRI r Q^s ( ryr,mRl\ , /■'^t 2\ ^orx 

To evaluate this expression, the renormalization matrix for the three-quark vertex T 
is needed in the mRI scheme. Therefore we perform a perturbative expansion of the 
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dimensionally regularized three-quark vertices: 

Tf- = rr + ^ rf^if^^Pk) + ^ ^ r^:r(/i,p;^)- (36) 

If we apply the projectors introduced in eq. (fT4l) and make use of their linearity, we 
find 

Comparing with eq. (|28l) reveals the identities 

= ^r^;T(/^,P.))' (38) 

which have to be evaluated at /x^ = (pf + P2 + Ps) / for the momentum geometries 
used in the simulations. Inserting this into eq. (|35l) yields the result for the scheme 
matching matrix of the vertex in first order: 

^ii W = + {vp, , ^rf™(M,Pfe)) + ^("J- (39) 

Together with eqs. (|34|) and (1331) this determines our scheme matching for the three- 
quark operator. 

The perturbative calculation of Tf^(n,pk) for the different three-quark operators is 
carried out in Euclidean space-time with off- shell quarks and gluons in Landau gauge 
^ = 0. It results in lengthy expressions and not all occurring integrals over the Feyn- 
man parameters can be solved analytically in closed form. The final evaluation of the 
integrals over the Feynman parameters and the construction of the projectors Pj as well 
as the evaluation of {vp. , ^r'''j5'(At,pfe)) were performed with Mathematica for all required 
momentum combinations. 

Note that it is important to exercise care when evaluating the terms of order 0(e°) in 
the three-quark vertex that define Tf^. Generally speaking one has to keep track of all 
possible terms in the Dirac structures proportional to e that multiply a 4 divergence, 
since these produce additional contributions to the scheme matching matrix. We have 
used the following strategy to treat the continuation to d dimensions. In a first step 
we have written our irreducible three-quark operators in four dimensions as linear 
combinations of the following basis operators: 

^pX.,t^.,u. =eciC2C3(^Ai ■ ■ ■^A,Ma)ci(C7p75)a/3(^Mi ' --D^^Up),, 
X {D^-^ . . . D^^dr)c3, 
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X {D^^...D^„{-f5d)r)cs, 

X {D,,...D,Sl^d)r),,. (40) 

As usual a, (3, 7 and r denote spinor indices, Aj, /ij, z/j as well as p and /x are space-time 
indices. Note that the operator W is equal to the operator V up to the position of the 
down quark. For the special case of operators without derivatives we have also used 

{C{-i)apf,)^pdpc2{l5u)rc3, (41) 
to access the operators of sub-leading twist. 

Then we have rewritten the three-quark vertices belonging to the above operator basis 
in terms of dimensionally regularized loop integrals. The corresponding Feynman di- 
agrams consist of three quark lines and one gluon exchange, leading to three strings 
of gamma matrices in the associated amplitudes. In the vertex two of these strings get 
contracted due to the presence of the {C . . .)a/3 structure. We can evaluate the remain- 
ing contractions of space-time indices using the li-dimensional Dirac algebra with an 
anticommuting 75 and the relation 

-C^,C-' = (7^)*. (42) 

This allows us to identify all contributions that are proportional to e°. Once these are 
determined we construct the regularized vertices of the irreducible three-quark opera- 
tors from the linear combinations of the A, V, W and U operators that were derived 
at the very beginning. Finally we evaluate the projections of eq. (|39l) resulting in the 
desired scheme matching matrices 2^r,Ms^mRi 



3.3 Renormalization Group Behavior 



Knowing the scaling behavior of Z provides a valuable consistency check of the re- 
sults. To this end one can compare lattice results that were derived at different renor- 
malization scales with the perturbatively expected scaling. This will be done in sub- 
section |4]2l 

Generally, the scaling behavior of a quantity is described by the renormalization group 
equation (RGE) and the related beta- and gamma-functions 

/3 = fi^-f^asi^i), (43) 
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^ = -Z-\^^)^^^—Z{^^). (44) 

Both 13 and the anomalous dimension 7 can be written as an expansion in the strong 
coupling: 



47r 



00 



Arc 



i+2 



i+1 



(45) 
(46) 



In the MS scheme, as well as the beta function are known to high order in perturba- 
tive QCD. The gamma function is operator dependent and in most cases not known to 
the same accuracy. 

Let us focus on the behavior of the renormalization matrix for the three-quark op- 
erators under a change of the renormalization scale. It follows from the RGB that a 
renormalization matrix can be converted from one scale fx to any other scale fl with the 
help of a scaling matrix AZ^^: 

Azf (^)zf (/i) = Azf (/i)zf (/i). (47) 

The matrix AZ^^{fi) again depends on the beta- function and the operator- specific 
gamma-function. We derive the scaling matrix for the three-quark operators by inde- 
pendently converting the renormalization matrices of the three-quark vertex and the 
three quark fields, cf. eq. (fTTl ). As the behavior of the quark fields is known to higher 
accuracy, we hope to get a better description for the three-quark operator renormal- 
ization by treating this contribution, just as the quark field scheme matching, to order 
al 

Hence we apply the expression (IA.4I) in Appendix IA.2I to convert the wave-function 
renormalization Zg together with its anomalous dimension from Appendix lA.3[ A lead- 
ing order formula is used for the scaling of the three-quark vertex renormalization 
matrix Z^'^ifi): 

AZf(^) = (a,(/i)-^«/2/^o)^^.. (48) 

We will use this to rescale our results obtained at a set of different fxs to the final scale 
fi = 2 GeV. 

We still need the gamma function 7^" of the three-quark vertex. Omitting contributions 
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that vanish for e — > 0, eq. (l36l) can be cast in the forni 

rr = rr + ^7.. (' - In ^) rr + + oial), m 

where is some momentum square, ji the scale of the as expansion and Cj is a finite 
term. We thus find for the renormalization matrix in the MS scheme 

,7r,MS/ \ r gnif^y ~ 2 , 4\ /;rn\ 

^ij W = Sij - -J^^iq + ^idR)- (50) 
Using eq. (f44)) and gndj,)'^ = g'^iJL~^ + 0{gj^ yields the anomalous dimension matrices: 

7S = -^7, + (^(«D- (51) 

The results for the anomalous dimensions are summarized in Appendix |Bl 

The whole procedure can only be expected to work in a "window", where the renor- 
malization scale is large enough for perturbation theory to be a good approximation 
and small enough so that the lattice cutoff is sufficiently far away. Hence we expect the 
condition 

Aqcd < /i' < — (52) 

di 

to restrict to a reasonable range. 



4 Lattice Calculation and Error Estimation 

4.1 Details of the Lattice Setup 

We work with non-perturbatively 0{a) improved clover- Wilson fermions and the pla- 
quette gauge action. We use gauge configurations generated by the QCDSF/UKQCD 
collaborations with two dynamical flavors of sea quarks, Uf = 2. The gauge is fixed to 
Landau gauge by minimizing the functional 

F^[U] =Y.Tr{u^{x) + U^ixy) =Y.RcTrU;:ix), (53) 

X,fl X,IJ, 

with 

U^(x) = G{x) Uf,{x) G{x + jiy. (54) 
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Table 1 

The used gauge configurations of the QCDSF/UKQCD collaborations. The first column shows 
the available lattice couplings /3, the second column the lattice volume and the following 
columns summarizes the hopping parameters k. 



/? 


V 


K 


5.20 


16^ X 32 


0.13420 


5.20 


16^ X 32 


0.13500 


5.20 


16^ X 32 


0.13550 


5.25 


16^ X 32 


0.13460 


5.25 


16^ X 32 


0.13520 


5.25 


24^ X 48 


0.13575 


5.25 


243 X 48 


0.13600 


5.29 


16^ X 32 


0.13400 


5.29 


16^ X 32 


0.13500 


J. 29 


z4 X 4o 


0.13jjO 


5.29 


243 X 48 


0.13590 


5.29 


243 X 48 


0.13620 


5.40 


243 X 48 


0.13500 


5.40 


24^ X 48 


0.13560 


5.40 


243 X 48 


0.13610 


5.40 


243 X 48 


0.13640 



Table 2 

Critical hopping parameters Kc and inverse lattice spacings in units of the Sommer parameter 
vq extrapolated to the chiral limit. 



p 


K/Q 


ro/a 


5.20 


0.13605 


5.454 


5.25 


0.136237 


5.880 


5.29 


0.136439 


6.201 


5.40 


0.136685 


6.946 



For any gauge link [/^(x) connecting the sites x and x + ft one can construct gauge 
equivalent links U^{x) by applying SU{3) transformations G{x). Landau gauge is re- 
alized by the gauge configuration U^{x) that minimizes the functional F'^[U] H16I17II . 
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We have performed computations on the ensembles of gauge field configurations listed 
in Table [U Table [2] provides further information on the chiral limit of these lattices. 
Our approach uses momentum sources [12J which improves the statistics significantly, 
as it "averages" over all space-time coordinates and thus is superior to point-source 
methods. Due to this benefit an order of ten gauge configurations per ensemble provide 
sufficiently high statistics for our purposes. Anyhow, the statistical error will turn out 
to be negligible compared to the systematic error related to the perturbative scheme 
matching and scale conversion. 

All calculations for the renormalization coefficients were performed on the Regensburg 
QCDOC machine with partitions of up to 128 nodes. The code implementation was 
done in C-i-i- using QDP-i-i-, the SciDAC data-parallel programming interface. We have 
taken over the lattice action and an even-odd-preconditioned stabilized bi-conjugate 
gradient method for the inversion of the Dirac-operator from the Chroma software 
library f,18J9J . 

For the discretized versions of the first and second covariant derivatives acting on a 
fermion field inside a three-quark operator we have used 

D^H^) = \ {U,{xmx + fi)- Ulix - fi)^ix - /i)) , (55) 

DlDli){x) = ^ {U^{x)U^{x + fi)ij{x + /i + z>) 

— Ufj_{x)Ul{x + fi ~ i>)'ip{x + fi — u) 

— Uji^{x — fj,)U^{x — il)ip{x — (1 + 0) 

+Ul{x - fi)Ul{x - fi - v)i){x - fi - v)) . (56) 



To minimize discretization effects, we choose the three external four-momenta pi of 
the quarks close to the diagonals of the Brillouin zone. Varying the squares of the 
momenta allows us to derive the renormalization coefficients at several different values 
of the renormalization scale /x. The Sommer parameter tq = 0.467 fm lfT4l is used to 
set the scale. 

At the end we have to perform a chiral extrapolation of the results. Therefore, we com- 
pute the amputated three-quark vertices T for fixed lattice size and fixed coupling P at 
typically three different values of the hopping-parameter k. For each k value we then 
determine the "renormalization matrix" k). The final mRI renormalization co- 

efficients are derived from a linear extrapolation in 1/k to the critical value 

1/kc. A typical chiral extrapolation is depicted in Fig. [T]and shows a reasonably linear 
behavior. 
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Fig. 1 . A typical chiral extrapolation for a renormalization coefficient of an operator with two 
covariant derivatives. Tiie curve shows a linear fit in 1/k at a fixed scale fi in the RI-MOM 
scheme, the cross marks the chiral limit. 
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4.2 Data Analysis 

The Jackknife method is used to estimate the statistical errors for all elements of the 
matrices /t). From the uncertainty in the At to the chiral limit we then 

determine the statistical error for each coefficient of the final result Z^^{ji). 

Due to its definition, may also contain imaginary parts before the scheme 

matching to MS. This can be seen best from eq. (|49|) . Applying the projector Pj to the 
right-hand-side of this equation yields the inverse of the mRI renormalization matrix. It 
is obvious that, apart from the real terms resulting from the expressions proportional to 
rf^^, we will also end up with a term proportional to PjCi from the finite contribution. 
Written in vector notation this projection reads {vp. , vc^ ) and this scalar product is not 
necessarily real. So finite contributions like Ci may render the mRI renormalization 
matrix complex. The scheme matching to MS should cancel these terms. However, 
since we perform a one-loop scheme matching only, we cannot expect a complete can- 
cellation of the imaginary parts. In practice it turns out that the residual imaginary 
parts are reasonably small and may be used as an independent estimate of the system- 
atic uncertainty induced by the scheme matching. Hence we will quote the real part of 
the renormalization matrices Z^^(n) as our results. 

The systematic error is estimated via the scaling behavior of the renormalization ma- 
trices. To this end we perform a lattice renormalization at ten different renormalization 
scales ji. We exemplify our procedure with Fig. [21 where we have plotted the scaling of 
one diagonal and one off-diagonal coefficient of the renormalization matrix for three- 
quark operators with two derivatives in the irreducible representation r^. The general 
behavior is typical for all representations: While the diagonal coefficient is of order 



17 



Fig. 2. Scheme matching and scale conversion for a typical diagonal (upper curves) and a 
typical off-diagonal (lower curves) renormalization coefficient Zij. Dashed curves: mRI(;u), 
dash-dotted: MS(Ai), solid: MS(2GeV). 




Fig. 3. Dependence on the renormalization scale for the same renormalization coefficients as 
in Fig. |2j however with different external momenta of the three-quaik operators. The meaning 
of the curves is the same as in Fig. [2l One finds good agreement between both momentum 
geometries. 
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one in all schemes, the mixing off-diagonal coefficient is close to zero. Let us now 
turn to the different curves in the figure. The dashed curve displays the results in the 
mRI scheme. As expected it shows a clear dependence on the renormalization scale: In 
the investigated range = 3 GeV^ — 100 GeV^ the renormalization coefficients vary 
by roughly 40%. The scheme matching to MS slightly shifts both mRI curves to lower 
absolute values (dash-dotted curve). In the third step we convert all results to the same 
renormalization scale fi = 2GeV (solid curve). This leads to a manifest flattening of 
the graph. Ideally the curve would be a constant now, but as already stated neither the 
lattice, nor the one-loop perturbative approach is free of systematic uncertainties in the 
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whole /i^ range. At low scales we expect the perturbative expansion to break down, as 
the coupling becomes too strong for our one-loop approach to hold. This behavior can 
be observed in the region ji^ < lOGeV^. On the other hand we also have to expect 
cutoff effects from the lattice calculation at large scales. Also this feature can be found 
in the figure. In the regime of fi^ = 10 GeV^ — 40 GeV^ the results appear almost flat, 
which may be interpreted as a scaling window. Note that although operators with a dif- 
ferent number of derivatives may come along with scaling windows shifted to slightly 
higher or lower values of yU^, all scaling windows seem to overlap in the quoted region. 
This indicates that the chosen approach for the non-perturbative renormalization and 
perturbative scheme matching is well justified. 

Each renormalization coefficient converted to the scale jl = 2 GeV can be reasonably 
described by the three-parameter fitting formula 

Zf = Ai, + Bij log a^/i^ + a,aV, (57) 

where ji stands for the original renormalization scale before the conversion. This for- 
mula incorporates 0{a^) effects as well as potential further logarithms in /i^ stemming 
from higher order corrections in the scheme matching and scale conversion. Our final 
result is read off at fi'^ = 20 GeV^. As our estimate of the systematic error we take the 
maximum of the differences with the values at 10 GeV^ and 40 GeV^. 

For the two renormalization coefficients plotted in Fig. |2]one finds in this way 

Zu = 1.3370(21)(252), Zgg = -0.0528(1)(171), (58) 

where the first error is the statistical and the second the systematic uncertainty. As 
expected for our one-loop calculation, the systematic uncertainty is larger than the 
statistical error and exceeds it by a typical factor of (9(10) to (9(100). It would be 
interesting to test our estimates for the systematic error by comparing them with future 
higher order calculations of the perturbative expressions. 

We have also investigated the dependence of the renormalization matrices on the ge- 
ometry of the external quark momenta p, . To this end we have calculated the Z matrices 
for modified angles between the three quark lines running into the vertex. In Fig. |3]we 
show a typical result for the same renormalization coefficients as in Fig.[2l We find that 
up to roughly 30 GeV^ even in the RI-MOM scheme no significant difference between 
both geometries occurs, whereas at still larger scales small deviations are observed. 
The scheme-matched curves in the MS scheme lie in most cases on top of each other 
and the final results are consistent within the systematic errors. For the renormalization 
coefficients derived with the modified momentum geometry, shown in Fig.[3l we find 

Zu = 1.3340(23)(146), Zgg = -0.0469(1)(250), (59) 

in good agreement with the values in eq. (l58l) . 
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5 Results for (2 GeV) 



Let us now present the renormalization matrices in the MS scheme at = 2 GeV for 
the different H(4) irreducible representations rf, r|, r-, rp-, of isospin-1/2 three- 
quark operators. For these representations we use the same notation as in Ref. [[TOll . 
The superscript denotes the dimension of the representation, with the underscore indi- 
cating its spinorial nature. The subscript distinguishes inequivalent representations of 
the same dimension. 

We will denote the statistical error by E^^ and the systematic error by E^^. As men- 
tioned, the latter is estimated by comparing the Z values at 20 GeV^ with those at 
10 GeV^ and 40 GeV^ Thus the final result reads 

Z,,±Ef^± Ell (60) 

As the statistical errors are almost two orders of magnitude smaller than the systematic 
errors, we will quote the statistical errors only in two explanatory examples and drop 
them for the rest of the paper. 

Another source of uncertainty results from the error of A^^. We use A*^^ = 261 MeV 
from Ref. [14J. Adding the two given errors we obtain a combined error of 43 MeV. 
While the resulting uncertainty is considerably smaller than E^^ for operators with zero 
and one derivative it becomes comparable with E^^ in the case of two derivatives. 

The mixing multiplets can be read off from Table [3} Multiplets of the same representa- 
tion and same dimension can mix under renormalization. Lower-dimensional operators 
can also mix into higher dimensional ones of the same representation by powers of the 
inverse lattice spacing i. We will summarize these ^ and ^ admixtures - wherever 
present - in the last columns of the related renormalization matrix (cf. the operator 
basis in Appendix [C])- 

We find that the results are similar for different values of the lattice coupling P and 
essentially identical at different lattice sizes, which again demonstrates the consistency 
of the approach. For reasons of better readability, we will only discuss the result for 
our finest lattice /5 = 5.40, x T = 24p x 48 in this section, while more details are 
given in Appendix O There we also give the operator bases explicitly. 

5.1 Representation rf- (Zero Derivatives) 

The following renormalization matrix belongs to the three-quark operators without 
derivatives in the irreducible representation rf . There are two multiplets that mix with 
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Table 3 

Irreducibly transforming multiplets of three-quark operators with isospin 1/2 sorted by their 
mass dimension. The subscripts /, g and h indicate that the covariant derivative(s) act on the 
first, second or third quark, respectively, cf. ifTOl . 





dimension 9/2 
(0 derivatives) 


dimension 1 1/2 
(1 derivative) 


dimension 13/2 
(2 derivatives) 


4 

-^r 


^(i),MA ^(i),MA 




.^(i),MA .^{j),MA MA 
//I ' //2 ' //3 ' 

.^(i),MA ^/^{i),MA ^/^{i),MA 
•-^g/tl ' *-^g?t2 ' ^gh3 


4 






.^(i),MA ^^(i),MA ^^(i),MA 
//4 ' //5 ' //6 ' 

.^(j),MA MA ..^(j),MA 






^^(j),MA 


.^(i),MA .^{j),MA ..^{j),MA 

\^ r pry • \^ p rci • p rrx • 

.^{j),MA .^{j),MA .^{j),MA 
^gh7 ' *-^g/i8 ' ^gh9 


12 


^W.MA 


.^(i),MA 
.^{j),MA .^(i),MA 


.^(i),MA ^^(i),MA ^^{j),MA ^^(i),MA 
^ffW ' '-^//ll ' '-^//12 ' ^ffl3 ' 

.^(i),MA .^(i),MA .^(i),MA .^(i),MA 
^ghlO ' ^ghll ' '-^g/il2 ' '-^g/il3 


12 




.^(i),MA ^^(j),MA 
,^(i),MA ^^(i),MA 


^^(j),MA .^(i),MA ,^(i),MA .^(i),MA 
' ^fn5 ' "-^//le ' ^ffl7 ' 

.^(i),MA .^(i),MA .^(i),MA .^(i),MA 
^ghlA ' '-^g/il5 ' '^ghW ' '"^g/il? 



each other and we find 



E 



sy 



0.6892 -0.0285 
-0.0065 0.6953 

1.7 X lO-'* 5.4 X lO-*^ 
2.6 X 10-' 1.5 X 10-4 

0.0151 0.0083 
0.0020 0.0163 



(61) 



The diagonal elements of Z are smaller than one and for both operators of almost equal 
size. The mixing off-diagonal elements are both negative and amount to roughly four 
and one per cent of the diagonal coefficients. Furthermore we see that the statistical 
error is of relative order 10"^ which renders it negligible compared to the systematic 
error. The latter one is two orders of magnitude larger and hence will be the dominating 
source of uncertainty in any renormalization of matrix elements on the lattice. 
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5.2 Representation r— (Zero Derivatives) 

There is only one multiplet of three-quark operators without derivatives belonging to 

12 

the irreducible representation r— . Therefore the renormalization matrix becomes one- 
dimensional. Our result reads: 



Again, the diagonal element is smaller than one and the statistical error is of the order 
10^^ so that the systematic error dominates. 

5.3 Representation r- (One Derivative) 

We carry on with three-quark operators with one covariant derivative. Since the statis- 
tical and systematic errors are of similar order of magnitude as for the above operators 
without derivatives, we do not quote them here anymore for reasons of better readabil- 
ity, but refer to the summary in Appendix O We obtain 



The diagonal component of the renormalization matrix for operators with one covariant 
derivative is larger than those for operators without derivatives. 

12 

5.4 Representation t— ( One Derivative) 

For leading-twist operators with one covariant derivative there are three multiplets 

12 

of the representation r— . One lower-dimensional multiplet that belongs to the same 
irreducible representation and is formed by operators without covariant derivatives can 
mix with them via one power of the inverse lattice-spacing ^, cf. Tableland ifTOl . We 




(62) 



Z 



1.1260 



get 



1.0540 0.1081 -0.0693 9.3 x lO'^ 



Z 



0.0564 0.9920 0.0483 -2.0 x 10"^ 



0.0033 -0.0028 1.0890 -2.1 x lO'^ 
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The coefficients in the last column, which describe the mixing with lower-dimensional 
operators, are rather small. 



5.5 Representation t§- ( One Derivative) 



This representation has four mixing multiplets. As for r- there is no mixing with lower- 
dimensional operators. Therefore these two representations are especially well suited 
for the evaluation of matrix elements involving three-quark operators with one deriva- 
tive. We find: 



1.0470 0.1066 -0.0675 -0.0013 
0.0544 0.9898 0.0487 -7.8 x 10"^ 
0.0080 -0.0064 1.0870 8.4 x lO'^ 
-6.4 X 10-^ -0.0026 0.0111 1.1320 



The mixing of the second multiplet with the first one is in the realm of ten per cent, 
while the mixing of the fourth multiplet can almost be neglected. This might be related 
to the fact that it contains three-quark operators with different chiralities of the quark 
fields than the other three multiplets. 



5.6 Representation Tf (Two Derivatives) 



We now proceed with leading-twist operators with two co variant derivatives. These 
operators belong to five inequivalent irreducible representations of the spinorial hyper- 
cubic group H(4). We start with the six multiplets of operators with two derivatives 
in Tf , which can mix with the two lower-dimensional multiplets, which contain three- 
quark operators without derivatives. The ^-admixture of the latter operators is seen in 
the last two columns of the renormalization matrix Z: 

[ 1.3390 0.0282 -0.0010 0.0306 -0.1620 0.0693 6.6 x 10"* 6.4 x 10"^ \ 

0.0167 1.2950 -0.0030 -0.0808 -0.0458 0.0750 -2.4 x lO""* -0.0040 

0.0022 0.0058 1.2710 -0.0019 -1.4 x IQ-"* 0.1167 7.4 x lO"' -0.0092 

0.0174 -0.0892 0.0468 1.3010 -0.0803 0.0464 8.7 x 10"^ 8.3 x 10"^ 

-0.0872 -0.0794 0.0708 -0.0564 1.2080 0.0615 4.0 x 10"* -0.0020 

\ 0.0249 0.0475 0.0550 0.0285 0.0618 1.2810 2.8 x 10"' -0.0063 / 

Here the statistical errors are by almost a factor of ten larger than for operators without 
derivatives. However, the relative systematic errors do not seem to change consid- 
erably. We furthermore observe a clear hierarchy with the diagonal elements further 
increasing relative to three-quark operators with one and without derivatives. 
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5.7 Representation T2 (Two Derivatives) 

This irreducible representation consists of six leading-twist multiplets. As by group- 
theoretical arguments mixing with lower-dimensional three-quark operators can be ex- 
cluded, r| is the representation of choice wherever possible when working with three- 
quark operators with two derivatives. The renormalization matrix at 2 GeV reads in the 
MS scheme: 

/ 1.2830 -0.0503 0.0142 0.0296 -0.2268 0.1382 \ 
-0.0287 1.2360 0.0137 -0.0994 -0.1214 0.1365 
-0.0065 -0.0132 1.2820 3.8 x 10"" -0.0089 0.1363 

Z = 

0.0329 -0.0138 0.0310 1.3690 -0.0637 -0.0021 
-0.0193 -0.0231 0.0253 -0.0679 1.3330 0.0412 
\ 0.0291 0.0207 0.0301 0.0035 0.0547 1.3770 / 



5.8 Representation r- (Two Derivatives) 

Again, six multiplets of operators with two derivatives contribute to this irreducible 
representation. They can mix with the lower-dimensional operators with one derivative 
of the representation r-: 

/ 1.3150 0.0220 -0.0133 0.0139 0.1879 -0.0742 0.0067 \ 

-0.0124 1.2560 0.0090 0.0882 -0.0621 0.0653 -0.0020 

-0.0084 -0.0176 1.2680 0.0017 -0.0123 0.1278 -6.6 x lO-"* 

Z = 

0.0349 0.0517 -0.0295 1.3240 0.0809 -0.0426 -0.0053 

0.0691 -0.0584 0.0565 0.0543 1.2020 0.0944 -1.6 x lO"'' 

-0.0410 0.0549 0.0507 -0.0291 0.0715 1.2890 -8.2 x 10"" / 

The mixing coefficients in the last column indicate that the ^ -contribution of the lower- 
dimensional multiplet might be important in practical applications. 



5. 9 Representation r— ( Two Derivatives ) 

Eight operator multiplets with two derivatives mix with the three multiplets with one 

12 

derivative and the operator multiplet without derivatives of r— . For reasons of bet- 
ter readability we split off the last four columns of the renormalization matrix, which 
describe the mixing with these lower-dimensional operators, and summarize the coef- 
ficients in a separate matrix Z'. Then the ^ mixing with the operators with one deriva- 
tives can be read off from the first three columns of Z', the ^ mixing from its last 
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column: 



z' 



1.3020 
-0.0212 
0.0112 
0.0013 
0.0749 
-0.0335 
0.0617 
-7.9 X 10-* 

-9.4 X 10-* 
-0.0059 
-0.0017 
0.0057 
0.0020 
-0.0044 
-0.0046 
0.0101 



-0.0515 0.0317 

1.2450 0.0210 

0.0113 1.2820 

-0.0046 0.0019 

-0.0050 0.0329 

-0.0137 0.0328 

0.0676 0.0472 

0.0089 -0.0043 



0.0039 0.0184 -0.1943 0.0940 

-4.3 X 10"" -0.1064 -0.1106 0.1174 

-5.7 X 10-4 0.0186 0.0579 0.0845 

1.3220 -0.0093 -0.0166 0.0071 

1.3470 -0.0019 -0.0359 

-0.0483 1.3000 0.0382 

0.0182 0.1031 1.2900 

0.0065 0.0027 -0.0024 



7.6 X 10"" 
1.0 X lO-"* 
0.0010 
0.0158 



-4.5 X 10-* 


0.0046 


0.0208 ' 


0.0060 


-0.0051 


0.0139 


0.0052 


-0.0033 


0.0025 


-0.0188 


0.0091 


-3.0 X 10'* 


-0.0041 


0.0070 


-0.0236 


0.0066 


-0.0084 


-0.0200 


0.0029 


-0.0075 


-0.0083 


-0.0046 


-1.8 X lO-'' 


-9.4 X lO-'' y 



-8.9 X 10 " \ 
7.9 X 10 " 
-4.4 X 10 " 

0.1027 
-3.2 X 10 " 
-9.3 X 10 " 
1.7 X 10" 
1.3560 I 



Compared to the diagonal elements around 1.3, the sixth operator mixes into the first 
with 0.2 and the ^ mixing contributes with a factor of up to 0.02. When renormalizing 
matrix elements with these coefficients, the mixing of the lower-dimensional lattice 
operators can only be neglected, if the matrix elements of these operators are consider- 
ably smaller than those containing operators with two derivatives. Whether this is the 
case or not must be decided in each individual case. As already stated earlier, whenever 
possible one should anyway make use of the operators of representation t|, because 
these do not suffer from mixing with lower-dimensional operators. 



5.10 Representation ( Two Derivatives ) 



This last representation also contains eight three-quark operators with two derivatives 
that can mix with the four lower-dimensional operators with one derivative of the same 
representation. For reasons of better readability we have again split off the last four 
columns of the renormalization matrix and display the related coefficients in the sepa- 
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matrix 


Z': 




















/ 


1.3340 


-0.0073 


0.0235 


-0.0019 


-0.0065 


-0.2075 


0.0934 


0.0061 


\ 






-0.0048 


1.2800 


-5.5 X 10- 


•* -1.1 X 10-* 


-0.0925 


-0.0853 


0.0902 


-3.2 X 10-* 








0.0124 


0.0037 


1.2850 


-1.2 X 10 -* 


-0.0060 


-0.0017 


0.1101 


-4.2 X 10-5 








0.0059 


0.0091 


-0.0024 


1.3290 


0.0014 


0.0078 


-0.0142 


0.1231 




z = 


























0.0301 


-0.0645 


0.0426 


8.1 X 10'^ 


1.3430 


-0.0877 


0.0288 


0.0013 








-0.0581 


-0.0559 


0.0530 


2.5 X 10 


-0.0606 


1.2850 


0.0649 


-0.0015 








0.0298 


0.0263 


0.0589 


-2.7 X 10-* 


-0.0027 


0.0605 


1.3380 


-3.3 X 10-* 






V 


0.0047 


-4.3 X 10"^ 


9.7 X 10 


0.0401 


-0.0015 


0.0043 


-0.0015 


1.3740 


/ 




/ 


-0.0017 


-0.0021 


5.5 X 10 '* 


-0.0026 \ 
















0.0116 


-0.0152 


0.0114 


1.1 X 10* 
















0.0025 


-4.2 X 10-^ 


0.0034 


-9.0 X 10-* 
















0.0028 


-1.1 X 10"* 


0.0068 


-7.7 X 10-5 












z' = 


























-0.0089 


0.0085 


-0.0066 


0.0040 
















0.0074 


-0.0057 


0.0037 


-1.0 X lO-"* 
















0.0083 


-0.0028 


0.0025 


0.0021 














I 


-0.0158 


0.0321 


0.0059 


-1.3 X 10-5 y 













Also in this case the renormalization matrix has mixing coefficients of a few per cent 
for the lower-dimensional operators. Whether or not they have to be taken into account 
depends again on the magnitude of the corresponding matrix elements. 



6 Renormalization of Moments of the Nucleon Distribution Amplitude 

Having presented the renormalization matrices in the previous section, the next step is 
their application to the moments of the nucleon distribution amplitude 

0'™"(/i)= Jldxi] x[x'^xl<P{xi,X2,x^,ii), (63) 
[dxi\ = dxidx2dxs6{l — Xi — X2 — x^). (64) 

For a detailed discussion of the calculation of these quantities in lattice QCD, the 
notation used and the physical interpretation of the results we refer to [8 9]. Here we 
focus on the behavior under renormalization and quote numbers only to exemplify the 
involved orders of magnitude. 

6. 1 Zeroth Moment 

The zeroth moment is linked to a matrix element of three-quark operators without 
derivatives that is proportional to /at, the normalization constant of the nucleon wave 
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function: 



{0\O%\p) = /^(zpi7i - iP2l2)N{p). (65) 

Here N denotes the nucleon spinor. The definition of O'^q can be found in [|9l. For 
the renormalization of the matrix element only the multiplet and the irreducible repre- 
sentation it belongs to is relevant. The first spinor component reads in our MA-isospin 
operator-basis: 

Ko), = -^^f'"' (66) 

12 

This operator belongs to the irreducible representation r— . Hence the bare value of Jn 
is renormalized by multiplication with Z{t^): 

= ^(rr) fN. (67) 
With fN/m% ^ 4.3 ■ lO^^ we find fW/^N ~ 3.5 ■ 10-^ cf. TableS 



6.2 Next-to-Leading Twist Constants Ai, A2 



For three-quark operators without derivatives we have also computed the renormaliza- 
tion matrix for next-to-leading twist. This enables us to renormalize the constants Ai 
and A2, which describe the coupling of the nucleon to two different interpolating fields 
used in QCD sum rules. Since we work in Euclidean space we have: 

(0|/:.(0)b) = -Aim^iV,(p), 
(0|A<,(0)|p) = -A2m^iV,(p). (68) 

We perform the following discussion in detail to clarify the treatment of mixing ma- 
trix elements and the required change of basis. In a first step we relate the operators 
sandwiched between the nucleon and vacuum to our MA-isospin basis: 

v3 

Hence Ai renormalizes like —8 O^"''^'^ and A2 like ^ O^i^'^^ . The operators O^'^^''^^ 
and O^'^'^^ form the basis for our renormalization matrix for the representation rf and 
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mix with each other under renormalization, i.e., introducing the abbreviations 

L2 = -lXi, (70) 

o 

we have 

L'f" = ^ij(rf) Li, 

Lr = Z2ii4)U (71) 

Substituting back the bare As yields finally the desired relation between the bare and 
renormalized values: 

A5^" = -Z(rf)2i A2 + ^(Tf)22 Ai, 

XT = Zirtl^ A2 - -| Z{Tt)u Ai. (72) 



(5. J First Moments 



The first moments of the nucleon distribution amplitude can be derived from matrix 
elements of three-quark operators with one covariant derivative: 

moXilp) = /^0'"^"[(Pi7i - P272)(*P373 - Eip)^,) - 2zpip27i72]iV(p). (73) 

The operators O^", where I + m + n = 1, are related to the MA operators of the 
representation r^— : 



),MAn 



{0Z), = l0T'- (74) 

The matrix elements of these operators yield the bare values for fN<P^^^, fN<P^^^ and 
fN<P^^^, as eq. ([731) shows. Using the definitions 



lOlO 



Ml = -() 
M2 = 

Ms = - 0^°°, (75) 
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Table 4 

Bare and renormalized values for (combinations of) the lowest moments of the nucleon distri- 
bution amplitude. In the case of the renormalized values the first error is statistical, while the 
second error estimates the systematic uncertainties due to renormalization and chiral extrapo- 
lation. The data are taken from 191 • 







Ai/GeV^ 


A2/GeV^ 


bare 


0.00429(7) 


-0.0729(14) 


0.1464(27) 


ren. 


0.00349(6)(12) 


-0.0498(9)(42) 


0.0985(19)(87) 




^100 


^010 


^001 


bare 


0.294(6) 


0.272(6) 


0.274(6) 


ren. 


0.346(8)(9) 


0.312(8)(13) 


0.314(8)(10) 




^200 


^020 


^002 


bare 


0.113(7) 


0.095(6) 


0.106(6) 


ren. 


0. 152(1 1)(83) 


0.127(10)(25) 


0.140(10)(17) 




^011 


^101 


^110 


bare 


0.065(4) 


0.069(6) 


0.071(4) 


ren. 


0.084(7)(31) 


0.112(9)(31) 


0.105(7)(2) 




Si 


S2 




bare 


0.840(10) 


0.722(16) 




ren. 


0.972(12)(13) 


1.021(28)(98) 






5100 


^010 


^001 


bare 


0.252(10) 


0.230(8) 


0.240(9) 


ren. 


0.370(14)(63) 


0.316(13)(16) 


0.336(14)(20) 



we obtain the following relation to the renormalized first moments of the nucleon 
distribution amplitude: 

/r 0°°''"" = Zir^^h. fM- (76) 
Upon using = Z{t^) In we get: 
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r''''' = -XjrZ{T,'"-^),M.. (77) 
Z{r-) 

Due to eqs. (|63l) and (t64l) the moments must comply with 

:= 0100 + 0010 + 0001 X 0000 ^ ^ (yg) 

While the bare values do not fulfill this equation (the sum equals 0.84), the sum of the 
renormalized values is found to be 0.97. Within errors this is in agreement with the 
constraint. 



6.4 Second Moments 



Analogously we can renormalize the second moments of the nucleon distribution am- 
plitude. A typical matrix element is given by 

+ zp3-E(p)7374(^Pi7i - ip2'l2)\N{p). (79) 

Again we investigate the multiplets and the representation to which the operators Ol^^, 
l + m + n = 2, belong by rewriting them in terms of our isospin mixed- antisymmetric 
basis. For the fourth spinor component we have, e.g.. 



(1),MA 
//6 



(1),MA 



(0-)^ = l^(9£r. (80) 



(1),MA 
gh6 
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With the new observables 



M[ = , 


^020 


M^ = , 


^002 


M^ = , 


^002 


Mi = . 




M^ = , 




M' =, 





^200^ 



iOll 



we find the renormalized quantities 



Here, the renormalized values must fulfill four constraints: 



5-100 






+ 




= 0100 


^010 


= /20 




+ 




= 0010 


^001 


= /02 




+ 




= 0001 




_ ^100 




+ 


^001 


= 1. 



(81) 



z{Tr) 

^020,ren ^ r^lDD^-. , w 

i002,ren ^ ry i DDArs , w 

= 12x [Z{^2 >5i-Z{T2 )6i)Mi, 

zijr) 

^ = ry, \2. Z[T2 )4iMi, 

0iio-- = -l^Z(rr-)5.M;. (82) 



(83) 

Comparing the results in Table |4] reveals good agreement within errors also for these 
quantities. This is an encouraging result, which demonstrates the consistency of the 
applied renormalization. It also makes us confident that we did not severely underesti- 
mate our errors. 



7 Summary and Conclusions 



In this paper we have set up a lattice renormalization scheme for three-quark operators 
based on the RI-MOM approach. The introduction of a suitable set of projectors fa- 
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cilitates the definition and practical evaluation of a renormalization matrix in our mRI 
scheme. In a second step we have performed a scheme matching to MS and a scale con- 
version to 2 GeV based on one-loop continuum perturbation theory. It was found that 
the statistical errors are negligible compared to the systematic uncertainties. Finally 
we have explained how to apply our results to the renormalization of low moments of 
the nucleon distribution amplitude. 
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A Perturbation Theory for the Three- Quark Operators 

A. 1 The Action 



In this appendix we summarize the conventions that we have adopted for the pertur- 
bative evaluation of the three-quark vertex V in one-loop order. We use the Euclidean 
action 

Se = j d'x[^{^,D, + m)^ + \f;,,F;, + ^^{d,Al){d,Al) 

-u''d^{d^5ab-gfab,Al)u% (A.l) 

with 

F;, = d,Al - d^Al + gfabcAlAl. (A.2) 

Here il) denotes the fermion fields, u the ghosts, the gluon field and Aq stands for 
the Gell-Mann matrices fulfilling [A^, Af,] = ifabc^c- The propagators in momentum 
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space are given by: 



pZ _|_ 



We work in the chiral limit m = and use the Landau gauge ^ = 0. 
A.2 The Scaling Function 

A three-loop formula for the function AZ^^{n) is given by: 



Az'-'i^r = o^si^r (1 + ^(71 - m + \ (^) ((71 - mo? 

+ 72 + ^i7o - Ml - ^270)) + 0{al), (A.4) 
with the abbreviations /3i = I3i/ (3q and 7, = 7j/(2/5o). 

A. i r/ze Quark Field Anomalous Dimension 



Finally we give the anomalous dimension of the quark field so that we can extract the 
scaling behavior from eq. (|A.4I) . It reads [|20ll : 

aj ra I 81^5 3 3 j 69 263 39 ^ 5 , 

+ (4^F'^-*'^''<"l4r"4'''^"8'-'^ 12"^ "32^ "l6^' 



Here 



C^ = 3, C'^ = ^, T = ^. (A.5) 
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B Anomalous Dimensions 



In the following we summarize the gauge invariant anomalous dimensions 7 for the 
three-quark operators. In the one-loop perturbative expansion they are identical to the 
anomalous dimension 7^", eq. (|5T| ). of the three-quark vertex in Landau gauge: 

7 = 7^k=0 + O(a')- (B.l) 



The operator basis is taken over from the linearly independent, isospin-1/2 H(4) ir- 
reducibly transforming multiplets of three-quark operators introduced in Appendix B 
of EQ] . We have checked that the eigenvalues of our anomalous dimension matrices 
agree with those presented by Peskin in ETI . 



B. 1 Operators without Derivatives 



We start with three-quark operators without covariant derivatives. For the irreducible 
representation rf we choose the operator basis 

To first order we find the anomalous dimension matrix 



(B.2) 




The operator basis for the representation rf is taken to be Oi = Oj . Then 



/ = S(2{-|). (B.3) 



B.l Operators with One Derivative 



In the following we summarize the anomalous dimensions for three-quark operators 
with one covariant derivative. The representation r- has the operator basis Oi = 

0){ With 

7^ = ^(2^-4). (B.4) 
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Z)12 

For the irreducible representation — we choose the operator basis 



Then the anomalous dimension is given by 



7 = 



47r 



16 
■ 9 



9 9 



2e-f y 



(B.5) 



In the basis 



Z)12 

of the representation T2 ~ we find the anomalous dimension matrix 



/ O £■ 22 



7 



47r 



2^- 9 

-I H- 






14 
9 



2e 



26 
9 



^ 



2^-4 y 



(B.6) 



B.3 Operators with Two Derivatives 



Finally, let us quote the results for operators with two covariant derivatives. In the 
H(4) -irreducible representation ^- the operators 



Ui — ^ffi , U2 — ^ff2 ) <^3 — L^//3 , 

n — n^i)MA ^ _ x.^(i),MA ^ _ ^^(i),MA 
- ^ghl : ^5 - ^gh2 ' ^6 - %/i3 ' 
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form a basis and have the anomalous dimension matrix 



7 



An 



1 

3 




5 
9 



2^-3 



4 
3 



8 
9 





10 

3 

5 
3 



16 

9 







2e 



41 



The same anomalous dimension matrix appears for the operators 



O4 



,^(i),MA 



(i),MA 
g/i4 ) 



O2 
O5 



,^(i),MA 



o 



(i),MA 
g/i5 ' 



O3 



,^(0,MA 



O 



(i),MA 
gh& ' 



belonging to t^^~. The representation r^^- with the operators 



04 



(i),MA 
//7 ' 
^/^(i),MA 



(i),MA 
//8 ' 
/-o(O.MA 



(i),MA 
//9 ' 
^^(i),MA 



has the following anomalous dimension: 



/o f _ 32 _2 



7 



47r 








9 3 
2^-3 





_ 10 
3 

5 
3 



16 
9 



2e-f 




O C 41 

2^- -9- / 



The irreducible representation r^^— and its operator basis 



O5 




O 



(i),MA 
//lO ' 
(i),MA 
ghlO ' 



O2 



,^(i),MA 
//II ' 
,^(i),MA 
*-^9/ill ' 



O3 
O7 



.^(0,MA 
//12 ' 
,^(i),MA 
*^3/il2 ' 



O4 



,^(i),MA 
^//13 
,^(i),MA 
*-^9/il3 
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has the one-loop anomalous dimension 



7 



4:71 



1 

3 





5 
9 






2^-3 - 
2^ 




31 
9 



2e 



5 
9 

















40 
9 



2e 



4 
3 





8 
9 






34 



10 

3 

5 
3 





16 
9 



29 



2^- 9 






41 











_4 

3 







2^ 



16 



(B.9) 



The same anomalous dimension can be shown to apply to the second twelve-dimen- 



sional representation r^f ^— and its operator basis 



O5 



,^{i),MA 
//14 ' 
,^(j),MA 



O2 



^/^(i),MA 



O3 
O7 



^^(i),MA 

,^(i),MA 
^ghl6 5 



O4 



^/^(i),MA 

/---)(«) -MA 
(-^g/ilT • 



C The Renormalization Matrices for Three- Quark Operators 



In this appendix we present the renormalization matrix Z{2 GeV) and the error ma- 
trices E^^. As the statistical errors are much smaller than the systematic uncertainties 
we do not quote them here. Note that the error due to the error of A'^^ leads to an 
additional non-negligible uncertainty (similar in size to E^^) for operators with two 
derivatives. We also list the operator bases for all irreducible representations. For the 
notation of the three-quark operators compare again [,10,| . The renormalized operators 
are related to the bare lattice operators by 

Op=%0,. (C.l) 



We have results not only for the two lattices presented in the following, but also for all 
other lattices in Table [TJ However, in order to keep the paper at a reasonable length we 
restrict ourselves to the largest lattices (24^ x 48) at f3 = 5.29 and (3 = 5.40. 
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C. 1 Operators without Derivatives in Representation Tf 



This irreducible representation contains two mixing multiplets. The renormalization 
matrix Zij is given in the following operator basis: 



We now present our chirally extrapolated results. 
/3 = 5.29, lattice size: 24^ x 48 



^ 0.6838 -0.0290^ 



-0.0066 0.6901 



0.0176 0.0095 
0.0021 0.0190 



(C.2) 



(C.3) 



(3 = 5.40, lattice size: 24^ x 48 



0.6892 -0.0285 
-0.0065 0.6953 

0.0151 0.0083 
0.0020 0.0163 



(C.4) 



12 

C.2 Operators without Derivatives in Representation t— 



Only one operator multiplet belongs to rp-. Hence, there is no mixing present and the 
operator basis is given by 



(i),MA 



(C.5) 



/3 = 5.29, lattice size: 24^ x 48 



Z = I 0.8047 ) , 
E'^ = ( 0.0176 



(C.6) 
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P = 5.40, lattice size: 24^ x 48 



Z = 1 0.8131 
= ( 0.0139 1 . 



(C.7) 



C.3 Operators with One Derivative in Representation r- 



In leading twist there is also no mixing for this representation. We take the basis 



(i),MA 
/I ■ 



13 = 5.29, lattice size: 24^ x 48 



13 = 5.40, lattice size: 24^ x 48 



(^1.108o) , 
( 0.0164 ) . 

Z^ (^1.126o) , 
E'^ = f 0.0172 1 . 



(C.8) 



(C.9) 



(C.IO) 



C.4 Operators with One Derivative in Representation Tf 
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There are four mixing multiplets and we take 



(C.ll) 



13 = 5.29, lattice size: 24^ x 48 

/ 



Z ^ 



1.0440 0.0989 -0.0653 9.1 x 10"^ 
0.0525 0.9841 0.0452 -6.6 x 10"^ 
0.0033 -0.0018 1.0730 -2.5 x lO'^ 
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'^0.0085 0.0190 0.0106 1.8 x 10"^^ 



0.0064 0.0127 0.0038 1.7 x 10"^ 
0.0013 0.0014 0.0136 1.6 x 10"^ 



P = 5.40, lattice size: 24^ x 48 

/ 



1.0540 0.1081 -0.0693 9.3 x 10"^ 
0.0564 0.9920 0.0483 -2.0 x 10"^ 
0.0033 -0.0028 1.0890 -2.1 x 10"^ 



^ 0.0062 0.0185 0.0095 1.5 x 10"^^ 



0.0061 0.0087 0.0045 1.5 x 10"^ 
0.0010 0.0017 0.0150 1.4 x 10"^ 



C.5 Operators with One Derivative in Representation T2 
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We work in the following operator basis: 

Oi = 0(;pA, o 



0.^0 



/5 ' 
(j),MA 
/7 ' 



2 
O4 



^/-,(i),MA 



o 



(i),MA 
/8 ■ 



(3 = 5.29, lattice size: 24^ x 48 
/ 



1.0350 


0.0975 


-0.0632 


-0.0011 


0.0502 


0.9813 


0.0456 


-6.7 X 10- 


0.0080 


-0.0065 


1.0720 


9.1 X 10-^ 


-4.8 X 10 


-0.0022 


0.0109 


1.1150 


0.0091 


0.0171 


0.0107 


0.0014 


0.0082 


0.0135 


0.0050 


7.0 X 10"^ 


0.0016 


8.0 X 10"^ 


0.0138 


4.9 X 10"^ 



1^8.3x10-^ 0.0015 6.9 x 10-^ 0.0174 j 



40 



P = 5.40, lattice size: 24^ x 48 



1.0470 0.1066 -0.0675 -0.0013 ' 
0.0544 0.9898 0.0487 -7.8 x 10"^ 
0.0080 -0.0064 1.0870 8.4 x 10"^ 
-6.4 x lO-'* -0.0026 0.0111 1.1320 ^ 

0.0093 0.0166 0.0096 0.0010 ^ 
0.0076 0.0094 0.0054 6.1 x 10"^ 
9.8 X 10-^ 6.4 X 10-^ 0.0146 3.8 x 10"^ 
\^5.8xl0 4 0.0013 8.0 x 10-4 0.0176 J 



(C.16) 



C.6 Operators with Two Derivatives in Representation rf 



Here, mixing with lower-dimensional operators occurs 



(i),MA 
gh2 1 



O2 - Off2 , 



0.^0 



(i),MA 
//2 ' 

(i),MA 
g/i3 ) 



0.^0 



(i),MA 



O7 



//3 
1 ^^(i),MA 



(i),MA 



(C.17) 



f3 = 5.29, lattice size: 24^ x 48 



/ 1.3280 
0.0142 
2.2 X lO-'* 



0.0296 
1.2900 
0.0055 



-0.0063 0.0281 -0.1538 0.0644 6.1 x lO " 7.3 x lO " \ 
-0.0070 -0.0772 -0.0416 0.0712 -3.3 x 10"" -0.0034 
1.2640 -0.0044 -0.0014 0.1108 5.9 x 10"^ -0.0080 



0.0228 -0.0754 0.0427 

-0.0809 -0.0729 0.0667 
0.0296 0.0536 0.0499 



1.2930 -0.0557 0.0339 0.0011 

-0.0473 1.2040 0.0536 4.0 x lO"" 
0.0330 0.0721 1.2560 3.1 x lO^^ 



V 

/ 0.0341 0.0165 0.0125 0.0163 0.0269 0.0101 3.2 x lO"'* 3.7 x lO"'* \ 

0.0210 0.0227 0.0075 0.0123 0.0265 0.0173 3.1 x 10-* 0.0025 

0.0072 0.0041 0.0148 0.0111 0.0055 0.0041 8.3 x lO"' 0.0057 

0.0203 0.0479 0.0123 0.0026 0.0733 0.0395 5.5 x 10 * 2.9 x 10 * 

0.0278 0.0273 0.0162 0.0238 0.0193 0.0249 4.7 x 10"^ 0.0033 

\ 0.0137 0.0137 0.0134 0.0123 0.0310 0.0469 1.3 x 10"* 0.0052 J 



7.2 x 10-" 

-0.0014 
-0.0053 J 



(C.18) 
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(3 = 5.40, lattice size: 24^ x 48 



z = 



1.3390 


0.0282 


-0.0010 


0.0306 


-0.1620 





0693 


6.6 X 


10-" 


6.4 X 10-" \ 


0.0167 


1.2950 


-0.0030 


-0.0808 


-0.0458 





0750 


-2.4 X 


10-" 


-0.0040 


0.0022 


0.0058 


1.2710 


-0.0019 


-1.4 X 10-* 





1167 


7.4 X 10-5 


-0.0092 


0.0174 


-0.0892 


0.0468 


1.3010 


-0.0803 





0464 


8.7 X 


10" 


8.3 X 10-" 


-0.0872 


-0.0794 


0.0708 


-0.0564 


1.2080 





0615 


4.0 X 


10-" 


-0.0020 


0.0249 


0.0475 


0.0550 


0.0285 


0.0618 


1 


2810 


2.8 X 


10-5 


-0.0063 / 



/ 0.0372 





0205 


0.0099 


0.0196 





0327 


0.0134 


3.6 X 


10-" 


4.3 X 10 " 


0.0228 





0179 


0.0070 


0.0123 





0303 


0.0195 


3.6 X 


10-" 


0.0023 


0.0071 





0052 


0.0098 


0.0108 





0058 


0.0031 


7.3 X 


10-5 


0.0055 


0.0208 





0493 


0.0138 


0.0033 





0732 


0.0390 


5.7 X 


10-" 


3.6 X 10 " 


0.0273 





0281 


0.0162 


0.0246 





0196 


0.0243 


5.9 X 


10-5 


0.0032 


V 0.0138 





0147 


0.0139 


0.0120 





0304 


0.0452 


7.7 X 


10-5 


0.0051 



C. 7 Operators with Two Derivatives in Representation T2 



This is the only irreducible representation of leading-twist operators with two deriva- 
tives that is not subject to mixing with lower-dimensional operators on the lattice: 



(i),MA 
gh4 1 



02^0 



(i),MA 
//5 ' 
(i),MA 
ghb 1 



(i),MA 
//6 ' 
(i),MA 



(C.20) 



(3 = 5.29, lattice size: 24^ x 48 



/ 1.2810 -0.0332 0.0047 
-0.0204 1.2450 0.0051 
-0.0080 -0.0125 1.2710 
0.0367 -0.0044 0.0288 

-0.0157 -0.0213 0.0235 



0.0284 
-0.0922 
-8.8 X 10-" 
1.3540 

-0.0610 



-0.0408 

1.3260 

\ 0.0322 0.0242 0.0269 0.0049 0.0622 

/ 0.0022 0.0467 0.0229 0.0165 0.0363 0.0418 \ 

0.0260 0.0379 0.0222 0.0071 0.0442 0.0476 

0.0050 0.0053 0.0142 0.0093 0.0195 0.0057 

0.0098 0.0280 0.0077 0.0243 0.0566 0.0362 

0.0079 0.0031 0.0041 0.0150 0.0170 0.0187 

y 0.0099 0.0100 0.0063 0.0042 0.0205 0.0683 J 



0.2071 0.1238 \ 
0.1045 0.1201 
0.0123 0.1309 
-0.0173 
0.0319 
1.3480 J 



(C.21) 
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(3 = 5.40, lattice size: 24^ x 48 



/ 1.2830 -0.0503 0.0142 0.0296 -0.2268 0.1382 \ 

-0.0287 1.2360 0.0137 -0.0994 -0.1214 0.1365 

-0.0065 -0.0132 1.2820 3.8 x 10"'' -0.0089 0.1363 

0.0329 -0.0138 0.0310 1.3690 -0.0637 -0.0021 

-0.0193 -0.0231 0.0253 -0.0679 1.3330 0.0412 

\ 0.0291 0.0207 0.0301 0.0035 0.0547 1.3770 / 

/ 0.0053 0.0525 0.0249 0.0118 0.0464 0.0477 \ 
0.0281 0.0416 0.0245 0.0101 0.0497 0.0508 
0.0040 0.0056 0.0153 0.0066 0.0145 0.0081 

0.0099 0.0257 0.0061 0.0283 0.0539 0.0366 
0.0070 0.0020 0.0029 0.0134 0.0176 0.0186 
\ 0.0093 0.0103 0.0063 0.0042 0.0209 0.0682 J 



(C.22) 



C.8 Operators with Two Derivatives in Representation r- 



This irreducible representation mixes with one lower-dimensional operator. The basis 
is 



Oi = 



(j),MA 
//7 ' 

(i),MA 
ghS 1 



02 = 



(i),MA 
//8 ' 

(j),MA 
gh9 ' 



o. = o 



(j),MA 
//9 ' 



07= - O 
a 



(i),MA 
/I ■ 



04 = 



(i),MA 



(C.23) 



13 = 5.29, lattice size: 24^ x 48 



/ 



1.3080 0.0156 



-0.0055 



-0.0113 1.2590 5.8 x 10"* 
-0.0069 -0.0177 1.2620 
0.0426 0.0415 -0.0296 
0.0635 -0.0557 0.0551 



0.0165 0.1784 -0.0734 0.0054 \ 

0.0785 -0.0553 0.0638 -0.0020 

0.0033 -0.0130 0.1216 -4.6 x 10"" 

1.3140 0.0541 -0.0238 -0.0056 

0.0487 1.2000 0.0841 -5.4 x 10"" 



\ -0.0471 0.0606 0.0468 -0.0322 0.0848 1.2620 -8.7 x 10"" / 



/ 


0.0190 


0.0300 


0.0345 


0.0119 





0095 


0.0099 


0.0039 


\ 




0.0085 


0.0300 


0.0328 


0.0204 





0150 


0.0106 


4.3 X 10 


4 




0.0060 


6.1 X 10-" 


0.0135 


0.0097 





0083 


0.0058 


5.1 X 10 


4 




0.0163 


0.0205 


0.0028 


0.0259 





0507 


0.0422 


0.0011 






0.0148 


0.0076 


0.0032 


0.0047 





0219 


0.0285 


0.0013 




v 


0.0101 


0.0086 


0.0113 


0.0050 





0248 


0.0489 


5.2 X 10 


V 



(C.24) 
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(3 = 5.40, lattice size: 24^ x 48 



/ 1.3150 0.0220 -0.0133 0.0139 0.1879 -0.0742 0.0067 \ 

-0.0124 1.2560 0.0090 0.0882 -0.0621 0.0653 -0.0020 

-0.0084 -0.0176 1.2680 0.0017 -0.0123 0.1278 -6.6 x 10"" 

0.0349 0.0517 -0.0295 1.3240 0.0809 -0.0426 -0.0053 

0.0691 -0.0584 0.0565 0.0543 1.2020 0.0944 -1.6 x lO"" 

y -0.0410 0.0549 0.0507 -0.0291 0.0715 1.2890 -8.2 x lO-"* / 

/ 0.0207 0.0371 0.0389 0.0154 0.0148 0.0080 0.0037 \ 

0.0073 0.0303 0.0387 0.0213 0.0194 0.0095 5.2 x 10-* 

0.0050 0.0026 0.0171 0.0076 0.0040 0.0077 3.8 x 10-* 

0.0141 0.0172 0.0042 0.0291 0.0447 0.0371 0.0016 

0.0132 0.0053 0.0027 0.0040 0.0150 0.0250 0.0015 

\ 0.0095 0.0087 0.0118 0.0053 0.0245 0.0485 6.4 x lO""* / 
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C.9 Operators with Two Derivatives in Representation t— 



Here, twelve multiplets of operators mix with each other under renormalization. Four 
of them have lower dimension: 



Oi = 


^fflO ' 


O2 


^^(i),MA 
~ ^//ll ' 


O3 


— ^ffl2 1 


O4 


— '-^//13 ' 


05 = 


^ghlO ' 


Oe 


_ /o(0,MA 


O7 


_ A^(i),MA 


Os 


_ (-nil) MA 
~ ^ghl3 ' 


09 = 


1 ^^(i),MA 

a ■' 


Oio 


_ 1 .^(i),MA 
— - *-^/3 ) 

a ■' 


On 


1 ^^(i),MA 

a ■' 


O12 


1 ^^(i),MA 
~ ~2 ^7 



We have split off the last four columns of the renormalization matrix, which describe 
the mixing with the lower-dimensional operators 09,...,Oi2, and display the related 
coefficients in a separate matrix Z'. 

p = 5.29, lattice size: 24^ x 48 



z = 



/ 


1.3020 


-0.0310 


0.0186 


0.0039 


0.0223 


-0.1658 


0.0771 


-0.0010 


\ 




-0.0162 


1.2540 


0.0113 


-4.3 X lO-'' 


-0.0941 


-0.0925 


0.1026 


8.0 X lO-'' 






0.0130 


0.0171 


1.2710 


-6.1 X lO-'* 


0.0206 


0.0630 


0.0720 


-3.3 X lO"'* 






3.6 X lO-'' 


-0.0057 


0.0024 


1.3150 


-0.0095 


-0.0177 


0.0080 


0.1001 






0.0771 


0.0010 


0.0308 


8.9 X 10-^ 


1.3350 


0.0166 


-0.0444 


-3.4 X lO-"* 






-0.0306 


-0.0114 


0.0310 


-1.7 X lO-'' 


-0.0403 


1.2920 


0.0306 


-9.1 X lO-"* 






0.0646 


0.0720 


0.0436 


6.6 X 10 "' 


0.0227 


0.1103 


1.2640 


2.0 X 10 "* 




\ 


-0.0010 


0.0085 


-0.0043 


0.0120 


0.0064 


0.0027 


-0.0021 


1.3330 


/ 
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0.0016 0.0116 0.0515 0.0352 5.1 x 10" 



2.2 X IQ-'* 0.0291 0.0459 0.0321 4.4 x 10 ' 
5.4 X 10-** 0.0041 0.0195 0.0249 3.2 x 10"' 



0.0216 0.0576 0.0399 

0.0138 0.0479 0.0329 

0.0077 0.0191 0.0182 

0.0015 0.0028 0.0014 0.0256 0.0018 0.0029 0.0026 0.0090 

0.0029 0.0086 0.0067 3.5 x lO"'' 0.0253 0.0378 0.0158 2.9 x 10"^ 

0.0066 0.0030 0.0029 5.4 x lO"'* 0.0152 0.0128 0.0177 4.9 x lO""* 

0.0045 0.0099 0.0104 9.8 x 10 ** 0.0146 0.0203 0.0511 4.4 x 10 ^ 

8.7 X lO"'' 0.0023 4.7 x 10 ** 0.0042 0.0019 0.0013 0.0021 0.0325 



/ -7.4 X lO-'' -5.6 X 10 ** 



Z' : 



-0.0057 
-0.0019 
0.0070 
0.0017 
-0.0049 
-0.0043 
0.0100 



0.0048 
0.0043 
-0.0186 
-0.0033 
0.0075 
0.0024 
-0.0041 



V 

/ 9.9 X lO-'' 6.5 X lO-'' 9.8 
0.0017 0.0055 



0.0051 
-0.0052 
-0.0028 
0.0088 
0.0068 
-0.0089 
-0.0075 
6.2 X lO-'' 



0.0189 \ 

0.0128 

0.0023 
-3.1 X 10"* 
-0.0249 
-0.0205 
-0.0088 
-0.0011 



5.4 X lO-'* 
0.0048 
0.0012 
0.0014 
0.0024 
y 0.0049 



0.0041 
0.0077 
0.0029 
0.0027 
0.0034 
0.0052 



X 10- 
0015 
0021 
0049 
0018 
0012 
0018 
0038 



0.0142 \ 
0.0097 

0.0024 
2.6 X 10"* 

0.0061 

0.0071 

0.0016 
4.3 X 10"* 



(C.27) 



(3 = 5.40, lattice size: 24^ x 48 



^ 1.3020 -0.0515 

-0.0212 1.2450 

0.0112 0.0113 

0.0013 -0.0046 

0.0749 -0.0050 

-0.0335 -0.0137 

0.0617 0.0676 

-7.9 X 10"** 0.0089 



0.0317 0.0039 
0.0210 -4.3 X lO"'* 



1.2820 


-5.7 X 10- 


0.0019 


1.3220 


0.0329 


7.6 X 10-** 


0.0328 


1.0 X lO-'* 


0.0472 


0.0010 


-0.0043 


0.0158 



0.0184 -0.1943 
-0.1064 -0.1106 

0.0186 0.0579 
-0.0093 -0.0166 

1.3470 -0.0019 
-0.0483 1.3000 

0.0182 0.1031 

0.0065 0.0027 



0.0940 -8.9 X 10"* ^ 

0.1174 7.9 X 10 ** 

0.0845 -4.4 X 10-'* 
0.0071 0.1027 

-0.0359 -3.2 X 10"* 

0.0382 -9.3 X 10"* 

1.2900 1.7 X 10 ** 
-0.0024 1.3560 j 
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/ 0.0192 0.0595 0.0395 9.2 x lO"'' 0.0111 

0.0132 0.0463 0.0328 5.4 x 10"^ 0.0274 

0.0078 0.0189 0.0197 2.2 x 10 0.0051 

0.0014 0.0030 0.0011 0.0167 0.0011 

0.0015 0.0065 0.0042 4.0 x lO"'* 0.0253 

0.0061 0.0014 0.0031 7.2 x lO''* 0.0146 

0.0044 0.0081 0.0096 0.0011 0.0130 

4.9 X lO "* 0.0021 5.5 x 10'"* 0.0061 



0.0533 

0.0433 
0.0195 
0.0031 
0.0322 

0.0146 
0.0170 
0.0015 8.7 X 10-"* 0.0014 



0.0372 5.4 X 10-* \ 
0.0315 1.7 X lO-'* 
0.0253 4.0 X 10 * 
0.0030 0.0077 
0.0126 1.8 X lO-'* 
0.0159 3.4 X lO-"* 
0.0481 3.4 X 10 ^ 
0.0310 



Z' 



( -9.4 X lO-'* -4.5 X 10-'' 0.0046 0.0208 

-0.0059 0.0060 -0.0051 0.0139 

-0.0017 0.0052 -0.0033 0.0025 

0.0057 -0.0188 0.0091 -3.0 x lO""* 

0.0020 -0.0041 0.0070 -0.0236 

-0.0044 0.0066 -0.0084 -0.0200 

-0.0046 0.0029 -0.0075 -0.0083 

y 0.0101 -0.0046 -1.8 X 10-* -9.4 x lO""* j 

X 0.0106 \ 

X 10-^ 0.0072 
0.0014 
1.4 X 10-* 
0.0061 
0.0038 
0.0020 



/ 0.0011 6.4 X 10 "' 9.3 
0.0013 0.0047 8.4 
6.3 X 10"* 0.0033 0. 



0.0049 
9.8 X 10 ^ 
0.0016 
0.0018 

y 0.0030 



0.0053 
0.0027 
0.0029 
0.0025 
0.0034 



0016 

0032 
0014 
0013 
0012 








0.0041 4.8 X 10-" j 
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C.IO Operators with Two Derivatives in Representation 



Finally, we have the representation T2 ~. As basis we take the operators 



Oi 


- ^fflA > 


O2 


_ .^(i),MA 

- (-^//is > 


O3 


_ cMi)MA 
- L^//16 ' 


O4 


_ ^^i:i)MA 
- ^ff l7 ? 




/-o(«).MA 
— '-^fffel4 > 


0, 


/-o(*).MA 
— '-^ff;il5 ' 


O7 


/'0(O.MA 


Os 


/'0(»).MA 
— ^ghl7 ' 


09 


1 ^^(i),MA 

" a ■ ' 






On 


1 .^(i),MA 

a ■' 


O12 


_ 1 ^(i),MA 

a 



Again we split off the last columns of our renormalization matrix, which describe 
the mixing with the lower-dimensional operators 0g,...,0i2, and display the related 
coefficients in a separate matrix Z'. 
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(3 = 5.29, lattice size: 24^ x 48 



/ 1.3240 
-0.0060 
0.0113 
0.0060 
0.0354 
-0.0526 
0.0314 
0.0051 

/ 0.0308 0. 
0.0117 0. 

0.0046 0, 
0.0013 0, 
0.0154 0. 
0.0211 0. 

0.0065 0, 

y 0.0015 0, 



V 



-0.0052 0.0194 

1.2740 -0.0029 

0.0057 1.2760 

0.0097 -0.0027 

-0.0538 0.0379 

-0.0480 0.0490 

0.0310 0.0553 

2.2 X 10-5 -8.2 X 10-5 0.0357 

0.0112 
0.0055 

0.0035 
0.0016 
0.0263 



-0.0020 
-4.8 X 10'5 
-7.2 X 10-5 

1.3210 
5.1 X 10'5 
3.0 X 10-" 
-2.0 X 10-^ 



0110 0.0057 0.0012 

0147 0.0026 1.4 x IQ-'' 

0050 0.0134 2.5 x 10"'* 

0025 0.0014 0.0201 

0369 0.0135 2.0 x IQ-'' 

0272 0.0138 3.7 x 10-'' 0.0170 

0126 0.0099 3.0 x IQ-* 0.0088 

0015 6.8 X lO-"* 0.0132 2.5 x 10" 



-0.0097 -0.2074 0.0980 0.0063 \ 

-0.0886 -0.0937 0.0962 -1.6 x 10'' 

-0.0060 6.5 X 10-5 0.1001 3.6 x 10-5 

0.0018 0.0073 -0.0138 0.1163 

1.3290 -0.0631 0.0136 0.0013 

-0.0532 1.2840 0.0509 -0.0014 

2.3 X 10-'' 0.0656 1.3130 -2.3 x 10-" 

-0.0015 0.0045 -0.0017 1.3470 

0.0406 0.0350 5.9 x IQ-" \ 

0.0530 0.0447 4.2 x 10-'' 

0.0091 0.0144 2.7 x 10"'' 

0.0029 0.0042 0.0131 

0.0694 0.0463 4.5 x IQ-'' 

0.0265 0.0450 6.7 x 10-" 

0.0180 0.0516 6.7 x 10"'' 

8.5 X 10-'' 4.5 X 10-'' 0.0423 



(C.30) 



Z' 



( -0.0021 -0.0019 

0.0118 -0.0143 

0.0028 0.0012 

0.0029 -0.0025 

-0.0092 0.0081 

0.0081 -0.0063 

0.0083 -0.0022 

Y -0.0178 0.0333 

/ 0.0031 0.0052 

0.0047 0.0081 
6.3 X 10 " 0.0036 

0.0032 0.0085 

0.0027 0.0043 

0.0015 0.0013 

0.0051 0.0045 

y 0.0041 0.0077 



4.7 X 10 " 
0.0118 
0.0027 
0.0093 
-0.0066 
0.0036 
0.0023 
0.0068 

8.6 X 10 " 
0.0021 
0.0030 
0.0065 
0.0025 
0.0011 
0.0013 
0.0055 



-0.0025 \ 
-1.4 X 10-" 
-2.8 X 10-" 
-7.5 X 10-5 

0.0041 
2.1 X lO-"^ 

0.0020 
-8.3 X 10-* J 

7.8 X 10"" \ 

0.0011 
7.2 X 10 " 
3.2 X 10-5 
2.1 X 10-" 
2.7 X 10"" 
7.4 X 10"" 
6.7 X 10-5 j 
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P = 5.40, lattice size: 24^ x 48 



/ 1.3340 
-0.0048 
0.0124 
0.0059 
0.0301 
-0.0581 
0.0298 
y 0.0047 



-0.0073 
1.2800 

0.0037 
0.0091 
-0.0645 
-0.0559 
0.0263 
-4.3 X lO-"* 



/ 0.0323 0.0107 0.0056 
0.0120 0.0169 0.0025 
0.0045 0.0053 0.0151 
0.0013 0.0028 0.0016 
0.0145 0.0347 0.0130 
0.0203 0.0274 0.0136 
0.0044 0.0111 0.0092 

y 0.0016 0.0013 0.0012 



0.0235 
-5.5 X 10-" 
1.2850 
-0.0024 
0.0426 
0.0530 
0.0589 

9.7 X 10-5 

8.8 X 10-" 
2.7 X 10-" 

5.9 X 10-5 
0.0124 

2.7 X 10-" 
2.4 X 10-" 
3.4 X 10-" 
0.0127 1 



-0.0019 -0.0065 -0.2075 0.0934 0.0061 \ 

-1.1 X 10-* -0.0925 -0.0853 0.0902 -3.2 x 10-" 

-1.2 X 10 " -0.0060 -0.0017 0.1101 -4.2 x 10 ^ 

1.3290 0.0014 0.0078 -0.0142 0.1231 

8.1 X 10-5 1.3430 -0.0877 0.0288 0.0013 



2.5 X 10 



1-4 



-0.0606 1.2850 0.0649 



-0.0015 



-2.7 X 10-" -0.0027 0.0605 1.3380 -3.3 x 10-" 
0.0401 -0.0015 0.0043 -0.0015 1.3740 



0.0103 


0.0416 


0.0375 


9.1 X 10"" ^ 


0.0025 


0.0533 


0.0461 


4.4 X 10"" 


0.0015 


0.0092 


0.0148 


2.0 X 10-" 


0.0017 


0.0018 


0.0030 


0.0126 


0.0279 


0.0628 


0.0431 


3.6 X 10-" 


0.0161 


0.0228 


0.0452 


4.2 X 10"" 


0.0079 


0.0133 


0.0487 


4.7 X 10"" 


7 X 10-" 


0.0013 


6.4 X 10-" 


0.0421 J 
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Z' 



-0.0057 
-0.0028 
0.0321 



( -0.0017 -0.0021 

0.0116 -0.0152 

0.0025 -4.2 X 10-5 

0.0028 -1.1 X 10-" 

-0.0089 0.0085 
0.0074 
0.0083 
y -0.0158 

/ 0.0018 0.0019 4.8 

0.0028 0.0060 0. 
7.4 X 10"" 0.0035 0. 

0.0016 0.0090 0. 

0.0014 0.0029 0. 

0.0020 0.0016 9.0 

0.0031 0.0032 0. 

y 0.0043 0.0043 0. 



5.5 X 10-" 
0.0114 
0.0034 
0.0068 
-0.0066 
0.0037 
0.0025 
0.0059 

X 10"" 6.3 

0011 8.7 
0027 7.2 
0065 2.6 
0016 1.4 
X 10-" 3.7 

0012 4.9 
0036 1.8 



-0.0026 \ 
1.1 X 10-" 
-9.0 X 10"* 
-7.7 X 10-5 

0.0040 
-1.0 X 10-" 

0.0021 
-1.3 X 10-5 

X 10-" \ 
X 10" 
X 10" 

X 10-5 

X 10" 
X 10" 
X 10-" 

X 10-5 J 



(C.33) 
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